Quasinormal modes for the scattering on a naked Reissner-Nordstrom singularity 
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What should be the quasinormal modes associated with a spacetime that contains a naked sin- 
gularity instead of a black hole? In the present work we address this problem by studying the 
scattering of scalar fields on a curved background described by a Reissner-Nordstrom spacetime 



with \q\ > m. We show that there is a qualitative difference between cases with 1 < q jt. 
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and cases with q r/m > 9/8. We discuss the necessary conditions for the well-posedness of the 
problem, and present results for the low damped modes in the low I and large I limit. We also 
consider the asymptotically highly damped quasinormal modes. We present strong evidence that 
such modes are absent in the case of a naked Reissner-Nordstrom singularity, corroborating recent 
conjectures relating them to classical and quantum properties of horizons. 



PACS numbers: 04.25.dc, 04.30.Nk, 04.70.Bv 



I. INTRODUCTION 

The naked Reissner-Nordstrom (R-N) singularity is a 
classical general relativistic solution in electrovacuum. 
The solution is expected to have a very limited mean- 
ing, due to the fact that such singularities cannot be cre- 
ated neither by a gravitational collapse, nor by dropping 
a charge into the black hole. (According to the weak 
cosmic censorship conjecture general naked singularities 
should be prohibited in general theory of relativity, al- 
though there are indications that by including quantum 
effects the violations of the conjecture could be consid- 
ered [![.) Moreover, a naked singularity created from 
some exotic initial data conditions should become quickly 
neutralized (classically, or via quantum pair production) . 
Some results also indicate that if one considers electro- 
gravitational perturbations the R-N naked singularity be- 
comes linearly unstable @- However it was discovered 
that the scalar field scattering problem on such a singular 
background can be still well denned [3H8| , since the waves 
remain regular at the origin. (However, the back reac- 
tion of the given scalar field configuration might still ex- 
cite some of the unstable electro-gravitational modes and 
this would eventually lead to a breakdown of the pertur- 
bation approach.) Despite the nice regularity property 
of the scattering problem, the spacetime is non-globally 
hyperbolic and the time evolution of the fields is not 
unique 0, Hoj . This means one has to specify an ad- 
ditional boundary condition at the singularity to obtain 
a fully unique time evolution. Another way of seeing the 
problem is through the language of operators: one can 
understand the spatial part of the wave operator as a 
positive symmetric operator acting on a L 2 Hilbert space, 
and then obtain the scalar field dynamics through a suit- 
able positive self-adjoint extension of such a symmetric 



operator 0, (One "preferred" way in which such a 
self-adjoint extension can be always realized is through 
the so called Friedrich's extension [3J, which will also be 
the case of this paper.) Anyway, after uniquely speci- 
fying the dynamics, one should be able to characterize 
the scattering by a set of characteristic oscillations, the 
quasi-normal modes. 

Low damped quasi-normal modes are in general used 
as a possible source of information about potential as- 
trophysical objects (such as neutron stars, black holes), 
and the highly damped modes are potentially interesting 
from the point of view of quantum gravity [TJ & 2] LJ Since 
a lot of work was devoted to the problem of quasinormal 
modes of the Reissner-Nordstrom black hole, it might be 
interesting to observe what happens if one transits from 
the R-N black hole case to the R-N naked singularity 
case (with a reflective boundary condition). Informa- 
tion about "what happens" shows how many features of 
the quasinormal modes of the black hole spacetimes are 
specific to the black holes themselves and what features 
survive much more general conditions. Let us also give 
one concrete example why asking what happens with the 
quasinormal modes of the naked R-N singularity might 
be interesting: In the black hole case the behavior of 
the asymptotically highly damped modes is widely sus- 
pected to be linked to the properties of the black hole 
spacetime horizon(s), (or more specifically it is consid- 
ered to carry information about quantum black holes). 
What does then happen with the asymptotically highly 
damped modes in case there are no horizons present? (If 
the behavior of the highly damped modes in case of R-N 
naked singularity would, for instance, resemble the be- 
havior of the highly damped modes of the black hole, it 
would be a disturbing fact from the point of view of the 
popular conjectures linking the modes to the horizon's 
quantum area spacing. Moreover, investigating the R-N 
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naked singularity from this point of view is especially at- 
tractive, as the R-N naked singularity is obtained by a 
continuous transition in the q, m parameters from the 
R-N black hole spacetime.) Thus, briefly, we hope that 
despite the fact that most likely the R-N naked singu- 
larity model does not correspond to a realistic physical 
situation, there are still many interesting things one can 
learn from such a model. 

The structure of this paper is as follows: In the second 
section we analyse the problem of the uniqueness of the 
time evolution of scalar fields on a R-N naked singularity 
background. In the sections three and four we analyse the 
properties of the effective potential for the scalar fields 
scattering and the geometrical optics (eikonal) limit of 
such a scattering problem. In the fifth section we de- 
fine analytically solvable potentials that can give good 
approximations to the problem of the low damped QNM 
frequencies (such that characterize the given scattering 
problem). In the sixth section we use those analytical 
approximations to derive semi-analytical results for the 
QNM frequencies in the eikonal limit. In the seventh 
section we use the numerical characteristic integration to 
obtain the low damped frequencies for the low values of 
I. In the eighth section we analyse what happens with 
the asymptotically highly damped modes and we suggest 
that in case of naked R-N singularity such modes do not 
exist (as one might exp ect considering some presently 
popular conjectures [ll|, Qjj])- We give the final conclu- 
sions in the section nine. We provide also appendices 
with more detailed results and some further suggestions 
for the analytical approximations of the problem. 



II. THE TIME EVOLUTION PROBLEM FOR A 
SCALAR FIELD IN THE R-N NAKED 
SINGULARITY 



In this section we will follow the standard analysis of 
the scalar field evolution in a curved background. (As 
an example of such an analysis see the treatment of 
Schwarzschild black hole perturbations in [li], EE] ■ For a 
review that presents also such techniques see for example 
16].) Take the Klein- Gordon equation for the complex 
(charged) scalar field: 



--d^v^d n*) = o, 



(1) 



with the metric line element given as 



g^dx^dx" = -f(r)dr + f{r)- L dr z + rW. (2) 

For the Reissner-Nordstrom (R-N) singularity the func- 
tion f(r) is in Planck units given as: 



1 2 

„. , 2m q 



(q 2 > m 2 ). 



(3) 



Take the decomposition of the field into the spherical 
harmonics 



^(t,r,6,(b)=J2Mt,r)Y T 



(4) 



After we separate the variables we obtain the following 
reduced equation 

d 2 Mt,r) f(r) d 



dt 2 



r 2 dr 
l(l + l)f(r) 



r Z f{r) 



dipi{t,r) 



dr 



(5) 



We are interested only in compactly supported data 
initial value problem: first assume that ipi (r, t) is ev- 
erywhere bounded and hence the Laplace transform of 
0z (r, t) exists: 



Jo 



(6) 



The Laplace transformed equation ([5]) gives the following 
equation: 



s 2 ^i(s,r) 



f(r) d 
r 2 dr 



7(r 



dTpi(s,r) 
dr 



l(l + l)f(r)~, . T . , 



where 



h(s,r) 



sipi(t,r) 



dipi(t,r) 



dt 



(7) 



(8) 



|t=o 



The solution that corresponds to the initial data term 
Ii(s,r) is obtained by the inverse Laplace transform of 
ipi(s,r) and the function ipi(s,r) is given as 



4>i{s,r) 



dr 1 Gi(s,r,r')Ii(s,r'). 



Here Gi(s, r, r') is a Green's function satisfying 



f(r) d 
r 2 dr 



r 2 f{r) 



dGi{s,r,r') 



(9) 



(10) 



2 , Ki + i)f(r) 



dr 



Gi(s,r,r , ) = S{r-r'). 



Since ij>i(r, t) was bounded, its Laplace transform ipi(s, r) 
must be also bounded in r. This translates to the bound- 
edness of the Green's function in r. The unique solu- 
tion of the given initial data problem is obtained if the 
condition of boundedness of the Green's function leads 
to a unique way to construct Green's function from the 
two linearly independent solutions Un, U12 of the homo- 
geneous equation 



fir) d 
r 2 dr 



r 2 /(r) 



dUn, 2 {s,r) 



l(l + l)f(r) 



dr 



Un, 2 (a,r) = 0. 



(11) 
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If f(r) goes to 1 at spatial infinity (the metric is asymp- 
totically flat), then there is only one solution of (|TT|) that 
stays bounded as r — >• oo. If there exists only one solution 
of the equation such that it is linearly independent 
from the solution bounded at infinity and in the same 
time it is bounded at r — > then these two solutions 
uniquely define the Green's function. If all the solutions 
of (|lip are singular at 0, there is no solution (bounded in 
r) of the given initial value problem such that it can be 
Laplace transformed (this can be taken as an indication 
that there is no solution at all). 

If both of the linearly independent solutions are reg- 
ular at 0, and (at least) two different Green functions 
lead to a function in the domain of the inverse Laplace 
transform, then there is no uniquely defined solution to 
the initial value problem. In the "worst" case there are 
infinitely many solutions, given by arbitrary linear com- 
bination of Un, U12, that are linearly independent to the 
solution bounded at r — > oo. In such a case the problem 
is underdetermined and one needs one more condition at 
r = that selects a unique Green's function between the 
different Green's functions marking different time evolu- 
tions. For each one of the choices of the Green's function, 
one can reproduce the calculation from [3, [HI and see 
that the quasinormal modes defined by 

1. the choice of the Green's function close to 0, 

2. the outgoing radiation condition, 

characterize the time evolution of the field at a fixed 
point within some specific time interval. Unfortunately 
for the case of R-N naked singularity (f(r) = 1 — 2m /r + 
q 2 /r 2 , q 2 > m 2 ) both of the linearly independent solu- 
tions Un, U12 are regular at and the problem is under- 
determined. It is easy to show that both of the solutions 
are regular at the origin. Write (jlll) as: 



2 d 2 C/ u , 2 (s,r) 
J[) dr 2 

2 , *0 + l)/(r) 



/(r)(2r-2m) 



dUn j2 (s,r) 



dr 



r* 



Un, 2 {s,r) = 0. 



(12) 




with 



"»-?±Vw + ^ eR - (15) 

hence both are regular. 

The fact that the problem is underdetermined is not 
surprising, since the space-time is not globally hyper- 
bolic and anything can fall out at any time from the 
singularity. This means the singularity has "hair" (car- 
ries some other information beyond the metric) and the 
quasi-normal modes obviously depend on the "hair" . Let 
us add here that, as we already mentioned in the intro- 
duction, the "hair" of the singularity relates to the ex- 
istence of many different self-adjoint extensions of the 
"Hamiltonian" operator in the equation. (In [9( one can 
find a nice analysis of the uniqueness of the self-adjoint 
extensions of such operators for many different types of 
naked spacetime singularities including the R-N naked 
singularity.) 

Is there any intuitive physical condition that we can 
further impose on the fields, that will uniquely select the 
appropriate Green's function? At least to get the ge- 
ometrical optics continuous extension of the black hole 
case one can impose the condition that nothing falls in 
or out of the singularity. This means there is neither ab- 
sorption nor superradiation in the scattering and the S- 
matrix of the K-G field is a unitary operator. What does 
this condition mean? The conserved current 4- vector for 
the complex Klein- Gordon field is given by 



J"(tf) = -igT (* V. - V, *) 



(16) 



Let us integrate the 4-current along a cyllindrical hyper- 
surface given by r = r and t G [ti, t 2 ] ■ Then since ((1 
is a conserved current it holds that 



Q(t 2 ) - Q(h 



dt d0 



(17) 



Now taking the r — > limit of the equation (|12j) one 
obtains the following 



d 2 Un,2{s,r) 
dr 2 



2m 



dUn i2 (s,r) 



1(1 + 1) 



dr 



Uu, 2 (s,r)=0. (13) 



This means the solutions Un >2 behave close to as 
Un, 2 (s,r) = exp(/3i, 2 r), 



Here Q(ti f2 ) is the integral along a hypersurface given 
by the interior of the cylinder at the constant time (ti or 
t 2 ), with a future orientated surface normal vector. E ro 
is a cylindrical hypersurface given by r = and \f—h is 



an induced density given as 
Furthermore is a normal 



r-Q sin(0). 
given as 

n (t,r,M) = (o, V7(r^),0,0). Since we want Q to remain 
constant with respect to time as we take the limit tq —> 
(nothing flows out or into the singularity), the equation 
(14) (TT71) reduces to the following 



lim 

ro->0 



dt d6 



n»J u 



-i lim /(r )r 



I 2 



dt 



d9 sin9 / d(f>[^d r y* - = 0. 



(18) 



4 



Let us analyse the equation (fT5)l in decomposition into 
spherical harmonics, hence 

$(f,r,M) = ^*(r,i)^(W). (19) 

I, in 

Then we can rewrite the equation (|18|) as 

lim /(r„)rg V / " d* [i/n (r, t)d r ^ (r, f) 
-#M)0rV>M)]| r=ro =O, (20) 

for arbitrary ii,i2. Since /(^o) r o <? 2 as r o 0, in 
order to fulfill the equation (|20[) we imposeH for every I 
and every t 

J r (*,0) - [^(r,t)W(r,t) - ipt(r,t)d r iPi(r,t)] lr=0 = 0. 

( 21 ) 

But this means that the functions ipi (r, i) should be al- 
ways constrained either by the condition ipi(0,t) = 
or d r ip(r, t)\ r= o = 0. (One might argue that it will be 
enough to claim that the fields and their first r deriva- 
tives should be real at zero, but that does not put any 
general constraint on the normal modes.) Now we obtain 
the function tpi (r, t) from the normal modes as 

/>oo 

^j(r,t) = / due- iwt [c 1 (u})Un(r,iuj) + C2(ui)Ui 2 (r,iu})], 
Jo 

and take the condition ^;(0,t) = for arbitrary time t. 
This translates to 




duj e- wt [ Cl (w)l7,i(0,«j) + CB(tt)l7B(0,tw)] = 0. 



(23) 

But since both £/ n (0,iw) = U a (0,iu) = 1 (see (HI)) ) 
and Fourier-like transform given by (1231) should not map 
non-zero functions to zero, one obtains the condition 
C2(w) = — ci(uj). This means we compose the relevant 
wave packet only from the following modes: 

Ui(r,iu) = U n (r,iu) - U l2 {r,iuj). (24) 

The same line of reasoning applies to the condition 
d r ip{r,t)\ r -Q = and the wave-packets that fulfill such 
condition must be composed entirely from modes given 
as (I > 0): 

- U n (iu},r) U l2 {iu,r) 

Ui(iu,r) = . (25) 

Pi P2 

For I = 0, we have the coefficient fa — and the solution 
Ui2(r, iuS) behaves as a constant for r f=s (with the value 



2 Such a condition on the radial part of the current at the origin 
(determining whether the scattering is absorptive, radiative, or 
superradiative) occurs in the formulation of [id ], where it gives 
constrains on the domain of the "Hamiltonian" providing sym- 
metricity of the operator. 



set to 1), which means that the condition (|2"5")) does not 
make strictly sense; in such case the modes are given 
simply by the function Ui2(r,iuj). 

This shows that the condition of S-matrix being a uni- 
tary operator gives additional constraints on the Green's 
function. The previous conditions mean that we shall 
consider only a linear space of wave-packets formed 
purely from modes that vanish at zero, or a linear space 
of wave-packets formed purely from modes whose first r 
derivatives vanish at zero. (Of course one cannot super- 
pose wave-packets formed from modes having the vanish- 
ing radial derivative at with wave packets formed from 
modes that are vanishing at 0, since the 4-current does 
not linearly depend on a wave-function.) 

So let us pick one of those two types of modes, fix the 
particular Green's function and employ the following rea- 
soning: Since Jo is conserved and the coefficients in the 
equation (|5|) are everywhere outside a smooth function 
of r, the function ij)(r,t), arising from compactly sup- 
ported initial data could be unbounded in time only in 
the case the wave-packet becomes concentrated around 
the singularity and slowly growing asymptotically with 
time into delta function 



lim i/)(r,t) ~ 6(r). (26) 



But this scenario is prevented by the boundary condition 
that "nothing flows in or out of the singularity at any 
time" . This means we expect the solution to be bounded 
with respect to both time and space and all the reason- 
ing based on the assumption of the existence of Laplace 
transform is justified. 

Further in the text we will employ the field vanishing 
condition at 0. (This boundary condition at the singu- 
larity corresponds to what is known as Friedrich's ex- 
tension of a symmetric operator.) Thus the quasinor- 
mal modes will relate to the scattering problem follow- 
ing from the time evolution determined by the boundary 
condition ip(0,t) = 0. But everything that we will do in 
the following text can be repeated for any other mean- 
ingful normal modes boundary condition giving another 
time evolution of the fields. The general dependence of 
the quasi-normal modes (characteristic oscillations) on 
such a boundary condition (given by some linear com- 
bination Un{r,iio) + K Ui2{r, iui)) describes the way the 
quasi-normal modes depend on the "hair" of the singu- 
larity. Since the quasi-normal modes are the ones that 
carry astrophysical information about the astrophysical 
sources, their dependence on the "hair" of the singular- 
ity (in some simplifying sense given by the complex pa- 
rameter K) has a potential astrophysical importance (of 
course in case naked singularities have any astrophysical 
importance) . 
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III. THE SCALAR WAVE SCATTERING ON A 
NAKED SINGULARITY 

Using fa defined as <f>i (r, t) = ripi (r, t) and x the tor- 
toise coordinate given by the condition: 



(27) 



one can rewrite the equation (|SJ) into the following form 

9 2 <MM) d 2 <f>i(x,t) 



dt 2 

with 

V(m,q,l,x) 



dx 2 



= V(m,q,l,x)<f>i(x,t), (28) 



1(1 + 1) 2m 



2q 2 



f(r(x)). 



r 2 (x) r 3 (x) r A (x) 
And, for the normal modes e~ tujt 4>i(r), we can write 
d 2 Mx) 



(29) 



dx 2 



+ [uj a - V(m, q, I, x)\ <f>i(x) = 0. 



(30) 



If \q\ > to , we can see that f(r) given by eq. (|3J) 
has no zeros for real arguments, but eq. (|27|) can still be 
integrated to give 



2m 2 — q 2 



arctan 



\J q 2 — m 2 
+ m ln(r 2 — 2mr + q 2 ) + C , 

where C is an integration constant. We remark here r« 
oo as r — ¥ co, but for r — > we have 

2m 2 — q 2 / —to \ 

a;(r — > 0) = — , arctan 



(31) 



\J q 2 — to 2 
TOln(<j 2 ) + 



\ \/ q 2 — to 2 / 
C = cc(0) = constant, (32) 



as can be seen in fig. ([T]). Further in the text we 
take the tortoise coordinate with the boundary condition 
x(0) = 0. It means that with respect to the usual (C = 0) 
tortoise coordinate such tortoise coordinate is shifted to 
the origin by the transformation x — > x(r) — x(0). 

The potential (j2U)) has for the ratio q 2 /to 2 less than ap- 
proximately 9/8 and the relevant x (in the naked singu- 
larity case the domain of x is constrained) 3 extrema, one 
smaller "outer" maximum, one dominant "inner" maxi- 
mum and minimum in the potential valley between them. 
(For r — > the function V(r) — > — oo.) For q 2 /m 2 more 
than approximately 9/8 the potential has only one max- 
imum (thus only one peak). These features of the poten- 
tial (I2U1) can be seen in the figure [U 

Moreover, in the case q 2 /m 2 less than approximately 
9/8 and / > 1, the inner maximum becomes completely 
dominant, making the outer peak negligible as compared 
to the size of the inner peak. For the outer peak and 
for I ^> 1 the terms 2™ — represent only a small 
correction (as compared to the term proportional to I 2 ) 
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FIG. 1. Typical examples of the behavior of the new tortoise 
coordinate defined by eq. (|31l) for a spacetime with m — 0.5 
and different values of q, taking C = 0. 



and the second peak will vanish in this approximation a10 
q 2 /m 2 =9/8. (This is because the outer peak lies always 
at r > 1. It is also quite obvious that the inner peak 

exists only due to the fact that the terms ^ — become 
dominant for r close to 0.) The features described in this 
paragraph can be observed in the figure [3J 

Further in the text we will use the following notation 
related to the potential parameters: By Vi,V2 we mean 
the heights of the two peaks. (Vi is the first larger peak, 
Va the second smaller peak.) By a% we call the curvature 
of the first peak and by a.2 the curvature of the second 
peak. Furthermore, by x\ max we mean the point of the 
location of the top of the first peak and by X2max the 
point of the location of the second peak. Many of these 
parameters can be exactly calculated in the I 3> 1 limit. 
The results of these calculations are given in the appendix 



IV. EIKONAL LIMIT - THE PARTICLE 
PICTURE 



First let us have a look at the geometrical optics 
(eikonal) limit of the Klein-Gordon equation, which is 
valid for I ^S> 1. In particular what does the geometrical 
optics limit tell us about the fundamental mode? The 
effective potential for the massless particle is 



V(r) 



f(r). 



(33) 



3 As trivially expected, we will see that the same result comes 
directly from the potential for the motion of a massless particle. 
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For 




9 

> - 



(36) 



there are no extrema. The maximum represents the un- 
stable orbit of a massless particle, and for the black hole 
case it is only the maximum that is relevant (because it 
is located above the horizon). For the naked singularity 
with q 2 /m 2 < 9/8 there exists also a stable orbit of a 
massless particle. However the fundamental mode in the 
geometrical optics limit is related to the unstable orbit 
as 



FIG. 2. Potential V(r) given by eq. ([29]) with I = 2, m = 0.5 
for q = 0.48,0.5,0.52,0.54 and 0.56. (The curves from left 
to right correspond to the increase of charge.) Note that the 
dashed part of the potential (for q — 0.48 and 0.5) is inside 
the black hole horizon. 



The effective potential goes to plus infinity as r — > and 
its extremes are given as 



Zm±rn\ 



ri,2 



(34) 



This means the potential has one minimum and one max- 
imum for 



m 2 < 8' 



(35) 



= \fV{r ma x) - 



V"{ 1 max )P{r m ax) 



2V(r r , 



(37) 



(For a very good paper that discusses this topic see [171].) 
Note that here \J ~ V ^ 7 ^^ 2 ^" 1 -^- gives the unstable 
orbit decay rate and furthermore holds the following 



V"( 1 max )/ (fmai) 



d 2 V(r) 
dr 2 



d 2 V(x) , 
dx 2 



(38) 

(Again x is the tortoise coordinate.) Particularly for the 
Reissner-Nordstrom black hole/naked singularity (with 
q 2 /m 2 < 9/8 ) equation (|37|l can be expressed as: 



2m 



18m 13g 2 + 24m 2 32mq 2 10q 4 



2r„ 



(39) 



.,, 3m+J$m 2 -8q 2 . , c 

with r max ~ , given as before. 

From the black hole QNM-s boundary conditions one 
can see that such massless particles decaying from their 
unstable orbit and eventually escaping either to the B-H 
horizon, or the infinity represent the fundamental QNM- 
s in the I ^ 1 (geometrical) limit. One can observe (see 
the discussion later) that the same holds in the naked sin- 
gularity case. In the I ^> 1 case the fundamental mode 
of the massless perturbations is described by the picture 
of a particle decaying from its unstable "photon" orbit. 
(Note that in the naked singularity case the stable pho- 
ton orbit makes the particle decaying inwards oscillate 
around the stable orbit.) This should suggest that in the 
geometrical optics limit 

1. one shall expect continuity of the fundamental 



mode (as a function of q, m) when turning from 
the black hole to the naked singularity case, 

2. for q 2 /m 2 > 9/8 there is a clear indication that 
there do not exist low damped quasi-normal modes. 

The way how to understand the continuity (for I 3> 1) 
of the fundamental mode as one transfers from the black 
hole case to the naked singularity case could be the fol- 
lowing: Consider the partial S- matrix Si(u>). Within the 
S-matrix one can look for poles on the complex plane: 
obviously any of the poles must be necessarily non-real. 
For the relativistic time dependence convention chosen as 
exp(— iujt) the poles with ujj > must be purely imagi- 
nary (as a result of the given Hamiltonian being a sym- 
metric operator) and they represent bound states with 
the non-relativistic energies given as E = —uSj. The poles 
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FIG. 3. Above: potential V(r) given by eq. (J29J) with m = 0.5 
for q — 0.52 (q 2 /m 2 < 9/8) and different values of I (same 
caption as on the right plot). Note in this case the existence 
of the secondary peak, reminiscent of the black hole potential. 
Below: the same, but this time for q — 0.6 (q 2 /m 2 > 9/8). 
There is no secondary peak in this case. 



with w/ < represent quasinormal frequencies. Some of 
the quasinormal frequencies correspond to the resonances 
in the phase factor related to the partial S-matrix, and 
the wave packets formed out of the resonance energies 
represent quantum particles tunnelling out of the poten- 
tial valley; for the non-relativistic quantum particles (for 
example) with the decay rate given as (w 2 )i = 2wflW/. 
Now in the naked singularity case the fundamental mode 
is represented by a resonance that has a low enough en- 
ergy (at the level of the smaller peak) and gets trapped 
in the potential valley in between the peaks for some 
time and then radiated away. The larger peak behind 
the smaller peak is hugely dominant (large I) and effec- 
tively acts to the wave packet as an infinite barrier. There 
are probably much higher resonance energies related to 
higher QNM frequency overtones that are determined by 
the higher peak, but left unaffected by the details of the 
small peak. So there are two different effective regimes 



for the quasinormal modes, in the first regime they are 
sensitive only to the details of the smaller peak, in the 
second regime they are sensitive only to the details of 
the larger peak. This splitting disappears in the case of 
small l-s since the two peaks are of comparable heights 
and effectively interfere. 

Let us finish this section with one more remark: In 
the black hole case one can relate the fundamental mode 
to the peak of the potential (independently of I ^> 1) 
following the way of thinking in [l8j]- The ratio of 
the amplitudes of the reflected and transmitted waves 
A re f(ui) / A tra ns{u) is for quasinormal modes 1. The same 
happens for real ui-s, when the energies ui 2 are close ("al- 
most at" ) to the peak of the potential. So for "almost" 
real w 2 (fundamental mode) the continuity of the ratio 
of reflected and transmitted amplitudes analytically con- 
tinued to the complex plane suggests that the real part 
of the frequency should be near the square root of the 
peak of the potential. (This intuition is then to some 
extend confirmed also semi- analytically, such as by the 
Ferrari and Mashhoon approach using the Poeschl- Teller 
potential [l9j . or by using the inverted harmonic oscil- 
lator potential.) In the case of a naked singularity life 
is not so easy (in the reduced ID problem we have no 
transmitted amplitudes), so this simple logic fails. (One 
can say luckily it fails as if it remained valid one would 
also for I 3> 1 expect a discontinuity in the fundamen- 
tal mode as one jumps from the black hole to the naked 
singularity. This is because the second, dominant peak 
suddenly appears in the domain of the x coordinate, due 
to the fact that the domain of x discontinuously jumps 
when passing from the black hole to the naked singular- 
ity case.) Certainly one can say that any of the quick 
fits, such as were done in the black hole case do not lead 
to anything close to the numerical data obtained in this 
paper. 



V. GENERAL ANALYTICAL METHODS TO 
MODEL THE PROBLEM 

In this section we want to suggest some ways to model 
the naked singularity scattering (with a general time evo- 
lution) by solutions of the equation with V(x, l 7 m, q) 
replaced by analytically solvable potentials. The reasons 
are the following: First, this section serves as a basis 
for the analysis of the I S> 1 cases provided in the 
section IVT1 Second, it demonstrates that the naked sin- 
gularity scattering is in principle treatable via analytical 
approximations and thus shows the general power of an- 
alytical techniques. (Also the approximations obtained 
here might produce some future results for the quasinor- 
mal frequencies. For some further calculations see the 
appendix [Bj Furthermore, for a very nice overview of 
the results for the quasi-normal frequencies and related 
transmission resonances of the analytically tractable po- 
tentials see [2lj|.) Third, it might bring more insights 
into the physics obtained through the exact scalar field 
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potential (|29|). 

The key point is to split the scalar field potential into 
different domains and approximate it on each of those do- 
mains (or directly the solution of the equation ([50)0 by a 
different, analytically solvable potential. Then one has to 
impose the standard procedure: The logarithmic deriva- 
tives of the solutions on different domains must be glued 
on the domain's boundary. We might add here that since 
the analytically solvable potentials are typically fitted by 
the parameters of the peaks one expects such approxima- 
tion to work for the modes that are not too damped (the 
low damped QNM modes). The approximation cannot 
be taken too seriously for the highly damped modes. 

Let us suggest modeling all the solutions of the relevant 
cases as follows: 

1. The case q 2 /m 2 > 9/8 for arbitrary I 

For this case it might be interesting to consider the 
Morse potential. (For the definition and origins of all 
the potentials used in this section see also [H|.) On the 
other hand the infinite valley behind the one and only 
peak may not be well modelled by only a finite valley 
that is behind the peak of the Morse potential. Thus 
we consider the solution in the region 'behind" Morse 
potential region to be given as (f2"4")l for x — > and the 
logarithmic derivative gluing condition originates from: 

1. < x < ax : use function given as (for the details 
see equation ([2"I|): 

4>(r) = r*(r) - Ar{e^ r - e?' r ). (40) 

2. a\ < x : use the appropriate solution of the equa- 
tion ([3"Uj) with the Morse potential: 

y ie -ai(x-x lma:c ) ^2 _ e -ai(a:-ii m „)j 

Here a\ is a point most conveniently chosen where the 
Morse potential is (for any arbitrary I there is such 
point for x > 0). This means a\ can be analytically 
given as 

ai = x lmax — ln(2). (41) 

Let us make here one remark: Instead of approximat- 
ing the original scalar field equation for r ss one can 
consider taking the equation (|30|) in the r — > approxi- 
mation. This leads to the following: 

The equation (g2]) has solutions <j> = Cix 1 / 5 + C 2 x 2/3 . 
The solution such that fulfills ^(0) =0 must have C\ = 
0. This approximation of the solution we think to be 
less exact as the approximation given as (1401) . as it is in 



fact effectively only the first non-zero term of the power 
series expansion of (|40p taken at 0. (For such reasons 
we decided to use in our analytical approximations the 
approximate solution (flUl) .) 

2. The case of q 2 /m 2 > 9/8 in the large I limit 

Let us consider now the limit of large I. In this case, 
since the peak already lies close to and the change in 
the highest power of inverted x becomes more and more 
dominant (with higher i-s), one might simplify the ap- 
proximation given in case IV ll by using the following po- 
tential: 

1. < x < x\ max : use function given as (|40[) 

2. Xi m ax < x: use the appropriate solution of the 
equation (150)) with the potential (Poeschl- Teller po- 
tential): 

V 1 cosh~ 2 [a 1 (x - x lmax )}. (43) 

Here we also used the Poeschl- Teller potential instead of 
the Morse potential, because the logarithmic derivative 
gluing condition is much easier to solve for the solutions 
of the Poeschl- Teller than for the solutions of the Morse 
potential. The expense of the simplicity is that the ap- 
proximation by such potential might be slightly less exact 
than by using the Morse potential, but one still expects 
it to be accurate enough. 

3. The case q 2 /m 2 < 9/8 for arbitrary / 

If one does not want to neglect the second smaller peak 
(so one is interested in resonances related to the valley 
between the two peaks), the analytically treatable poten- 
tial describing all the resonances (and also higher damped 
QNM frequencies) could be: 

1. < x < ay. use function given as (|4"U[) . 

2. a\ < x < a 2 : solutions of the equation (|30[) with 
the Morse potential 

3. a 2 < x: solution of the equation (|30[) with the 
Poeschl- Teller potential: 

V 2 COsh.~ 2 [a 2 (x ~ X 2 max)}- (45) 

Here a\ is best chosen as in the case IV II and a 2 is cho- 
sen to be such that the resulting potential is continu- 
ous. Generally a 2 has to be obtained numerically, but for 
larger i-s (certainly Z = 20 is more than enough, as we 
checked), the curvature of the second peak is very small 
comparing to the curvature of the first peak and also the 
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curvature of the first peak becomes (for large l-s) very 
large comparing to the scales of the potential, so that 
one can calculate (for I = 20 with a good approximation 
at least to the 6 decimal places) a 2 just by 

1. taking the Poeschl- Teller potential to be constant 
and given by the hight of the outer peak 

2. simplifying the Morse potential by the following ap- 
proximation: 

V(x) « 2V 1 e- ai{x ~ Xl — >. (46) 

The resulting formula then becomes 

a-2 = xi max — In ( — \ J . (47) 

ai \2ViJ 

Then one can easily show that a 2 — > as I —¥ 00. 



the valley between the two peaks becomes flat as com- 
pared to the difference between the height of the first 
peak and the value of the potential at the bottom of the 
valley. If the resonances related to the smaller peak lo- 
cate close to the top of the peak (they should since the 
first peak grows with I an makes the tunnelling harder), 
then one might effectively approximate the case IV 31 by 
the Poeschl- Teller potential: 

V 2 cosh' 2 [a 2 (x - x 2ma x)}, (48) 

for < x. (One can also derive this approximation 
straight from the logarithmic derivative gluing condition 
in the appendix IB 21 by taking the I — > 00 limit.) 

VI. THE LARGE I LIMIT - THE ANALYTIC 
APPROACH 

In this section we want to analytically confirm the re- 
sults obtained directly from the eikonal limit. 



4. The case of q 2 /m 2 < 9/8 in the large / limit 

Now consider what happens with the resonances re- 
lated to the valley between the two peaks in the large I 
limit. In such case the dominant peak grows to infinity 
as compared to the smaller peak, the infinite valley be- 
hind the dominant peak shrinks to 0. Also a\ — ¥ and 



1. Calculations for q 2 /m 2 < 9/8 

The case IV 41 from the previous section is easily solv- 
able. The solution in the region x > must fulfil the 
outgoing radiation condition. (To be exact after ana- 
lytically extending the solution to the complex plane it 
must give asymptotically the outgoing waves on the line 
(u)x)i = 0.) Such solution is given as follows 



<I>r{x) = C 2 e iux F 2 i (gu92,g 3 , (1 + exp(2a 2 (x - a; 2max ))) x ) , (49) 
where F21 is the standard hypergeometric function (of the type 2-1) and 



1 




v 2 _ 


= 2 + 




1 — 

ai 



9i = i; + \h-^ = -92 + l (50) 



93 = 1- — . (51) 
a 2 



The boundary condition at gives </>_r(0) = 0, leading to 



, 1 ix/V? 1 i\/V 2 iui , , , . 

F 2X - + -^A - - 1 , (1 + exp(-2a 2 x 2ma;c ))- i = (52) 

2 a 2 2 a 2 a 2 



Unless q 2 /m 2 is not too close to the upper limit given for I 
large as 9/8, one can reliably approximate e~ 2a2X2mal « 
and the condition ([52^1 turns to be: 



F 2 \ 



WV2 1 
a 2 '2 



ry 2 



a 2 



OL 2 



(53) 



Now (|53p can be rewritten through Gamma functions 
and reduces to the simple problem of finding poles of the 



product of Gamma functions 



r 



ioj — i\fV 2 
a 2 



r 



iuj + i\JV 2 
a 2 



0. (54) 



The poles are located at — n for n being a natural number 
and this gives 



u; = ± y \ 2 - HI 2 ! n + - 



(55) 
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This is precisely the formula for the lowest damped black 
hole QNM frequencies. This means in the I 1 limit we 
see (at least for the lowest modes) a continuous transi- 
tion from the black hole case to the naked singularity 
CclSG, clS expected. (At least in the case where it holds 
that e -2 " 22:2maa! s» 0, which becomes a less accurate ap- 
proximation for q 2 /m 2 close to 9/8. On the other hand 
for such ratios of q 2 / m 2 the peak is almost vanished and 
one can assume that the decaying circular orbit eikonal 
picture has already quite limited sense.) 

All this means that in the regime I 1 the lowest 
QNM frequencies should be located close to the first peak 
of the potential and are (up to certain n) given by the 
formula (|55[) . They represent resonances such that they 
are insensitive to the details of the larger peak as in their 
case the larger peak can be already seen to act effectively 
as an infinite potential barrier. Hence we can conclude 



Here Q is defined as Q — uj/(2a\). For not highly damped 
frequencies in the I ^> 1 limit it must hold that a) e IR. 
But by plotting the right side of the equation (|55|) minus 
the left side of the equation (|58|) in Mathematica it can 
be shown that the equation ([55)) does not have any real 
solutions. For more than large intervals of values of Co 
the absolute value of the left hand side of the equation 
minus the right hand side of the equation always seems to 
grow (almost) linearly with respect to Co from the value 
approximately given as 0.35 at Co — to infinity. It might 
not be hard to prove also analytically that (|55|) does not 
have real solutions by using many relatively simple prop- 
erties of Gamma functions. This means there are no low 
damped modes in the case of q 2 /m 2 > 9/8 and I ^> 1, 
as expected from the eikonal limit / particle picture. (It 
means the damping of the fundamental mode grows to 
infinity as I — > oo and grows no less rapidly than a%, 
hence no less rapidly than cubically with I.) 



that in the I ^> 1 case we are able to match the geomet- 
rical (eikonal) limit of the original problem considered. 

2. Calculations for q 2 /m 2 > 9/8 

Take the case lV2l from the previous section. On the left 
hand side of the logarithmic derivative gluing condition 
we need a logarithmic derivative of the function: 

r* L = C L r (e Pir - e &r ) . (56) 

On the right hand side of the logarithmic derivative glu- 
ing condition we need the logarithmic derivative of the 
solution on the interval [x\ ma x, oo]. It is given as in (|49l) . 
only the parameters are related to the first and only peak 
in this case. Then the logarithmic derivative gluing con- 
dition leads to the following: 



(58) 



I 

VII. THE NAKED SINGULARITY FOR THE 
SMALL WAVE MODE NUMBERS - NUMERICAL 
RESULTS FOR THE FREQUENCIES 

In [22j], the quasi- normal modes of a scalar field in an 
electrically charged Vaidya background were studied us- 
ing a similar numerical setup to what we are going to use 
here. An interesting investigation done in [12] tried to de- 
termine what would happen with the quasi-normal modes 
as the time-dependent background approached a naked 
singularity, but the numerical code used was not suitable 
for following the field evolution after the extremal con- 
figuration (q — m) was reached. In our present work, 
we have a static configuration that describes a naked a 
singularity in order to study the properties of the quasi- 
normal modes. 

Our objective in this section is to solve eq. (|28[) with 
potential (|29p numerically, in the case where q > m as de- 
scribed in the last section. To do this, we rewrite eq. (|2"5)l 
in terms of the light-cone variables u = t—x and v = t+x, 
where x corresponds to the tortoise coordinate (l27l) . as 

d 2 <f> d 2 ^ d 2 4> 

Wd^ = - 4 d^ = v{r) ^ (59) 



Vi 



r (- iu - 1 /1 Vi 1 A r ( %uJ 1 1 /I _ Vi 1 I 

i \ 2«i 2 Y I ^ ~r i) 1 V 2ai 2 Y 4 a* T 4 

If one takes the I — > oo limit then the equation (|57|) becomes 

r ^ + i7S + !) r (- M 



i 

2V6 



1 

2V6 



lUJ - 



1 

2V6 



2a~i 



he 



foe 



e P2 



— = + coth(^3) 
v3 
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that can be integrated with the boundary conditions 

cj)(r = 0,t) = <f>{u,v = u + 2x ) = 0, (60) 

(v-va) 2 

(f>(u = 0,v) = e 57^- , (61) 

where condition (|60|) is a necessary condition on the field 
<j> near the origin (see the discussion on fig. [5] below) 
and condition (f6"Tj) defines an "arbitrary" relevant initial 
signal to be propagated. We use the algorithm 



>N 



f»s 5 KA„A„, (62) 



see where A„ and A„ are the integration steps in u and u, 
respectively, and the definitions of cW etc can be seen on 
fig. [U Note that here V is the potential (f2T))) evaluated 
at the same r coordinate as tfis (and aW). 



r = j» o 
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FIG. 4. Example of the numerical grid used for integrating 
eq. (I59[) with boundary conditions (|60|l and (|61[) . The points 
marked with "x" are out of our domain, values of <f> at the 
positions marked with filled circles are given by the boundary 
conditions and the values at the empty circles are obtained 
with the algorithm (|62p . 

As we can see in fig. [SJ the boundary conditions (|6"tJ)) 
and (|61[) ensure the necessary conditions on the fields cj) 
and i\) near the center. As we discussed previously in 
section HH the physically correct boundary condition for 
ip is ip(0,t) = 0. From this we must have for 4>(r,t) — 
rtp(r,t) that 0(0, t) = and <j/(0,t) = 0. In fact, the 
quantity 



£ = 



(hp 
di 



d.r 



^7 + h£ +V(r)<fr 



dx 



(63) 



is invariant along the i-evolution governed by ([5Uf with 
the boundary conditions <p(r = 0,t) — <p(r — > oo,t) = 0. 
Taking into account ([2T]) , the integrand can be expressed 
near r = as 



^4 W 



(64) 



with the prime denoting the derivative with respect to r. 
Any initial condition with finite £ , as those ones obeying 
dSB and (ED, is such that 



r 2+a 



(65) 



for r — > 0, with a > 0, implying that <p'(r = 0,t) = 
for initial condition with finite E obeying <fi(r = 0, t) = 0. 
Moreover, we have checked that our boundary conditions 
for (j> also reproduce the <f>'(Q,t) = condition, by testing 
the code with the condition <f)(f, i) = 0, where r is a very 
small value, and we found no differences in the late time 
evolution of the fields. 
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FIG. 5. Above: Behavior of <f> with I = 2 as a function of r 
near the center r = for a late time = 350, shown here 
in order to exemplify the effect of conditions (|60[) and (|61[) in 
the numerical integration, for a spacetime with q = 0.5 and 
q = 0.52. Below: The same as in the left plot, but this time 
for the function ip = cj>/r. 

In the left plot of fig. [HI we present some typical time 
evolutions of <p, for a I — 2 and different q/m ratios. In 
the right plot we present the obtained frequencies of the 
QNMs (fundamental mode) in the ujr x ujj plane. We 
can see a discontinuity in the frequencies as q/m — > 1, as 
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was expected from the discussion of the potential V(r) 
(see fig. H]). We also point here that we see no significant 
changes, but rather a smooth behavior as q 2 /m 2 — > 9/8 
(q/m — > 1.06 in the plot). But we see a point of inflec- 
tion in u)r at q/m « 1.16, for which we did not find an 
analytical explanation. 

We remark here some numerical issues that prevented 
us from extending the results shown in fig. [5] for lower 
and higher q/m rations. As q/m — > 1 and <fi becomes 
less damped, different modes subsist for longer times and 
and longer evolutions are needed in order to obtain the 
clear frequency of the fundamental mode. On the other 
hand, as we increase q/m, <\> is damped so quickly that 
we cannot observe enough oscillation cycles to obtain the 
frequencies. These are the issues that have limited our 
results. 




oji changes significantly with q/m. In the upper plots 
(q 2 /m 2 < 9/8), \uji\ decreases with I, that is, the damp- 
ing time is longer. In the lower plots (q 2 /m 2 > 9/8), we 
have the opposite tendency. This behaviour is connected 
to the potential V(r) shown in fig. [3l It might be also 
interesting to mention that in case q 2 /m 2 > 9/8, there is 
a qualitative similarity in the behaviour of the imaginary 
part of the frequencies as a function of I, between the 
case when I is small and the I large limit. 
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FIG. 7. Above: frequencies of the fundamental mode as a 
function of I for m = 0.5 and q = 0.52 {q 2 /m 2 < 9/8). Be- 
low: same as above, but this time for m = 0.5 and q = 0.6 
(q 2 /m 2 > 9/8). 
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FIG. 6. Above: <f>(xF,t) with I = 2 at xf = 100 for a space- 
time with m = 0.5 and different values of q > m. Below: Fre- 
quencies of the fundamental mode with / = 2 in the lur x uji 
plane, parametrized by the q/m ration. 

Finally, in fig. [7] we explore how the frequencies of 
the QNMs change with /. As usual in black hole scat- 
tering problems, we see that the oscillation frequency 
lur increases with I. But the qualitative behaviour of 



VIII. ASYMPTOTICALLY HIGHLY DAMPED 
MODES FOR THE SCALAR FIELD IN THE R-N 
NAKED SINGULARITY SPACETIME 

Let us explore the following limit, I will be fixed and 
n — > oo, assuming that there exist an infinite number 
of QNMs. Assume further that w 2 , uji diverge as n — > 
oo. In such a limit, to have a good understanding of 
what the QNM boundary conditions are, it is convenient 
to take an analytic continuation of the solution into the 
complex x plane and impose the outgoing wave boundary 
conditions on the (Stokes) line Im(ujx) = 0. This is for 
almost imaginary asymptotic frequencies (uj ~ iuji — > 
—zoo asn-> oo) the line Re(x) = 0. The key point is to 
realize that, in the asymptotic case, uo 2 3> |^(i)| holds 
everywhere apart of a tiny region around 0, where the 
potential V(x) can be already approximated by 



V(x) 



2 

9x2 



(66) 



(There is one subtlety: The potential, and also later the 
solutions, are multi- valued functions in the complex vari- 
able, but after specifying a branch cut together with the 
branch, this does not cause any problems to our analysis.) 
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So let us make the following statement: The solutions of 
the equation 



d 2 Mx) 
dx 2 



2 
9X 2 



<M^) =0 



(67) 



are everywhere good approximations to the solutions of 
the scalar field equation for the asymptotically highly 
damped modes. (The region where the approximations 
x « r 3 and V(x) ~ l/r(:r) 6 cease to hold is the region 
where we can already neglect the potential as a whole 
with respect to the w 2 term.) 

The equation (|6"7) has a general solution given as 



ne equation (|fc> Y j) has a general solution give: 
(f>l(x) = ^Jujx [AJ 1/6 (u>x) + BJ_i/q{ux)] , 
re J„ is the Bessel function. The solution fi 



(68) 



where J a is the Bessel function. The solution for which 
tpi(x) fulfils the vanishing boundary condition at is 
given by B = 0. Our QNM boundary condition says that 
this solution taken along the Stokes line Im(uix) = 0, 
should give purely outgoing radiation for lux — > oo. But 
in the ujx — > oo limit we know that we can approximate 
the solution (via Bessel functions approximations in such 
a limit) as 



Ay/cox Ji/e(ux) « A-y^cos (uix — — ^ 



(69) 



.4 



2tt 



This linear combination is independent on lj and does 
not give the purely outgoing radiation (~ e lux ), which 
suggests that our QNM boundary conditions cannot be 
fulfilled for asymptotically highly damped modes. 

This can be taken as an evidence that the asymptoti- 
cally highly damped modes (loj — > — oo) do not exist for 
the R-N naked singularity and this can be taken as a con- 
firmation that the highly damped QNMs could link to the 
black hole spacetime horizon's properties (non-existence 
of the horizon leads to non-existence of the modes in the 
asymptotic limit). One can see also some link to the 
fact that for the extremal R-N black hole the spacing 
between the asymptotic frequencies goes to zero as pro- 
portional to the surface gravity (this might make a whole 
infinite tower of arbitrarily highly damped modes even- 
tually "collapse" to a single mode). Moreover, we assume 
that the non-existence of the highly damped modes can 
be proven along the same lines also for another naked 
singularity spacetimes (like negative mass Schwarzschild 
spacetime etc |9|). 



IX. CONCLUSIONS 

In this paper we analysed the problem of the scalar 
field scattering on a R-N naked singularity background 
from the point of view of quasi-normal modes. The evo- 
lution on the R-N naked singularity is non-unique un- 
less one specifies an additional boundary condition rep- 
resenting a "hair" of the singularity. The quasi-normal 



modes then carry information about the "hair". We 
applied a particular boundary condition, that nothing 
comes out, or in from the singularity and analysed ana- 
lytically, as well as numerically, the characteristic oscilla- 
tions of the scalar field perturbations (low damped quasi- 
normal modes). We analysed the eikonal I 1 case via 
the analytical approach confirming the intuition obtained 
through the massless particle viewpoint, and showed that 
an approach based on analytical approximations can be 
useful also for the small I wave mode numbers. For the 
small l-s we calculated the frequencies numerically via the 
characteristic integration method. We also suggested ar- 
guments showing that the asymptotically highly damped 
modes (limit I fixed and n — > oo) do not exist in case of 
R-N naked singularity. This might confirm the intuition 
one has about such modes from the black hole physics 
(and from the presently popular conjectures [ill. ll2jV 

The basic results can be summarized as follows: for 
the low modes and the large I there is a continuous tran- 
sition in the low damped QNM modes between the R-N 
black hole and the R-N naked singularity. However, when 
the ratio q 2 /m 2 becomes larger than approximately 9/8 
then the picture becomes significantly different and the 
low damped modes do not exist for large l-s. (This is a 
very different picture from the BH based intuition.) For 
the small / numbers the modes face a discontinuous tran- 
sition when transiting from the black hole to the naked 
singularity. Furthermore, the / dependence \loi\ (for small 
I) changes as q 2 /m 2 becomes larger than approximately 
9/8: \uji\ decreases for q 2 /m 2 < 9/8 and increases for 
q 2 /m 2 > 9/8. It might be interesting to notice that for 
q 2 /m 2 > 9/8 the increase of as a function of I (for 
small l-s) matches the behaviour of \uj]\ for large ^-s. In 



the case of large l-s and q 2 /i 



i 2 > 



9/8 we have shown 



that \loi\ of the fundamental mode grows at least cu- 
bically with I and thus, as we already mentioned, the 
low damped modes do not exist. For the asymptotically 
highly damped modes our results seem to suggest that 
they do not exist, which means that the imaginary parts 
of the frequencies are bounded. 
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Appendix A: Calculations of some of the important 
quantities 

In this appendix section we provide a list of some of 
the quantities (also some of their derivations) that occur 
in the calculations relevant to the problem analysed in 
this paper. Write V(r) in a convenient form: 



V(r) 



A 



B 



C_ 



D E 



(Al) 
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with 



1. The large I limit 



A = 1(1 + 1), (A2) 

B = 2m[l - 1(1 + 1)], (A3) 

C = g 2 Z(Z + 1) - (2m) 2 - 2q 2 , (A4) 

D = 6mq 2 , (A5) 

£ = -2q 4 . (A6) 

Now let us calculate the second derivative of the potential 
with respect to the tortoise coordinate: 



d 2 V(x) _ d 2 V(r) (fry + dV(r) d 2 r 



dx 2 



dr dx 2 



dr 2 \ dx 

At r max the derivative dV(r)/dr is zero, then necessarily 



d 2 V(x) 
dx 2 



d 2 V(r) 
dr 2 



(A8) 



J \x„ 



and since 
d 2 V(r) 



dr 2 

we obtain 
d 2 V(r) 



6A 12B 20C 30L> 42£ , k . 
— + — + — + — + — ( A9 ) 



dx 2 



2m q 2 



6A 12B 20C 30L> 42E 

T "T" 5 I a H~ « ~\~ o 



(A10) 



For the global extremum condition holds: 
dF(r) 



dr 



2Ar„, 



-^Br^ ax +4C^ max +5Dr max +6E = 0. 

(All) 

This equation can have maximally 4 roots. From (|2T))) 
it can be easily seen that for r — > oo (I2TH) becomes for any 
l,q,m positive (dominant term is 1(1 + l)/r 2 ) and goes 
to 0. In case of naked singularity for r — > the potential 
goes always to — oo (dominant term is —2q 4 /r 6 ), in case 
of black hole for r —> r + it always goes to (dominant 
term is f(r)). This means, together with (lAlljl . that 
for any m, q, I (black hole or not) there are always either 
two local maxima and one local minimum, or one local 
maximum without local minima. 



For I ^> 1 the following approximations hold: For large 
I the maximum of the inner peak can be approximately 
found analytically and behaves as: 



V3k 



vf+i) 

Also holds the following: For i<1 

,3 



(A12) 



r 



(A13) 

The parameters ot\, V\ turn in the large / limit to be: 
2 [1(1 + I)} 3 / 2 



(A14) 



W + i)? /2 

3 3 / 2 



9 



(A15) 



This implies also the following results: For I — > oo 





2 


(A16) 


Oil 


i 

~ 2^3' 


(A17) 


V-2 


~ 1(1 + 1), 


(A18) 


Vi 


->0, 


(A19) 


Oil 


->o, 


(A20) 




— > oo, 


(A21) 


V(x) 


2 r 

9ar 


(A22) 



Appendix B: Some suggestions for the analytical 
treatment of the scattering for small wave mode 
numbers 

Write the two linearly independent solutions on the do- 
main i in the form convenient for the logarithmic deriva- 
tive gluing condition: d(^n + K^m). If the domain 
extends to the infinity and is taken to be the asymp- 
totically incoming wave solution and the asymptoti- 
cally outgoing wave solution, then Ki(u) is the S- matrix 
(ratio of the coefficients of the outgoing and incoming 
waves). Anyway, the logarithmic derivative gluing con- 
dition at the boundary of the regions i and i — 1 turns to 
be: 



15 



Ki = 



<5 



2*(t-X)2 



VE "( t -l)2^2 



(Bl) 



This is the condition connecting the if, coefficient with 
coefficient through the values of the solutions and 
their derivatives at the boundary of the regions. On one 
side consider asymptotically free-wave region, whereas 
on the other side consider infinite potential barrier-like 
boundary condition for ^S. In such case one can through 
equations of the type (|B1[) connect (after finite number 
of steps) the S-matrix of the external (asymptotic) re- 
gion with the boundary condition at the origin (where 
we can conveniently put K\ to be zero). Then one gets 
the S-matrix expressed through a function of the bound- 
ary values (at all the boundaries) of the solutions and 
their derivatives. This will be algebraically more diffi- 
cult condition than the one coming just from posing the 
outgoing wave condition in the external region, as we 
did in the previous sections (case of large I and q 2 /m 2 
larger than approximately 9/8). The reason for posing 
the more complicated condition of the type (IB1[) is that it 
might be easier (and also brings more insight) to look for 



a real energy resonances (by plotting phase change and 
looking for rapid phase shifts, hence rapid change of the 
logarithm of the S-matrix), than to try to numerically 
solve the more simple outgoing wave condition (with the 
help of appropriate computer software) . Such resonances 
should then 1 to 1 correspond to the low damped quasi- 
normal frequencies. Let us apply these ideas and try to 
calculate the expressions for the partial S-matrix. 



1. The case for q 2 /m 2 > 9/8 

This is the case IV 11 We are gluing only two regions. 
The solution on the left side of the boundary l , such 
that fulfils *l(0) = 0, is again given as ([iD|. On the 
right side (the Morse potential side) the solution is given 
as: 



^ R (x) = Cx[* R/I {x) + S(u, a Xj Vi)* R/0 (x)] = Cie-4 [e"^, s 2 , z) + S(u, a u V l )e iux M(s' 1 , s' 2 , z)] . (B2) 

By ty R /o(x), ^r/j(x) we mean asymptotically outgoing and incoming wave solutions in the sense that ^r/j(x — ¥ 
oo) — > e^ tux and ^r/q{x — > oo) — > e lljjx exactly. Furthermore M(., ., .) is being the Kummer (confluent) hypergeo- 
metric function and 



1 



si = x -i- 

l Oil 



V\ iu) 



2iuj 
Oil ' 



S 2 = 1 



2iu) 

Oil 

ilJVi 



4iuj 



a i 



Oil 



exp(-ai(a; - xi max )). 



Also Si(oj, Oil, Vi) is the S-matrix. The logarithmic derivatives gluing condition at ai — x\ max — ln(2) is then 



1 

a i 



filial _ ^g^ax 



o 0iai 



u\z{ai) . s[ M (s[ + l,s' 2 + l,z(ai)) 
~^~ +tU+ 7 2 M( S ' 1;S ' 2 ,z( ai )) ' 



Considering that the following holds: 



(B3) 
(B4) 
(B5) 

i 

(B6) 



one can rewrite (IB6I) into: 



z(ai) 



i^y/Vi 

Oil 



(B7) 



f(ai) 



1 

a i 



i2y/Vi +iuj+ -j 



,m(V 1 + 1,4 + 1,^ 



M 



( ' / ,/ ^WI 

^*1'*2> ai 



(B8) 



The partial S-matrix can be obtained from (jBll) [Ki = 0) Since for real uj holds that ^r/o = ^*r/i an d is a 
as: real valued function, one can immediately observe that 











*'l 


R/O 
*R/0 
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(for lj G R) S*(uj) — S[ (u>), hence partial S-matrix is unitary. 



From (|B9[) we can conclude the following 



5/(w, V\,ax) = -e 



/(«i) 



/(«! 



J_ 4. Sl^ 1 



.02 »i 



J_ , /3 1 e3i°i-fee< 3 2°l 



A/ 



!l+1 , S2+1 ,2i^ 



s 2 . , / i±*JV\ 

Ml s 1 ,s 2 ,—^- — 



M [ St .si, , 



(BIO) 



2. The two peak case, q 2 /m 2 < 9/8 



The two peak case is a more complicated problem. Anyway, one can write the logarithmic derivative gluing condition 



as 



a 1 z(a 2 ) . Si(ui,ai, V )M{si,s 2 ,z(a 2 )} + e r 

VjJ- 



»M{«i,« / 2 ,«(a 2 )} 



S^uj, a u V )M{si,s 2 ,z(a 2 )} - e i2uJa2 M{s[, s' 2 , z(a 2 )} 
Si(u, ai, V Q )M{ Sl + 1, s 2 + 1, z(a 2 )}^ - e i2ua ^M{s[ + 1, s' 2 + l,z(a 2 )}^ 
S 1 {u,a 1 ,V )M{s 1 ,s 2 ,z(a 2 )}- e i2 ^M{s[, s' 2 , z(a 2 )} 



(Bll) 



lid 



F 2 i{gi,g 2 ,g 3 - 1, [1 + exp(2a 2 (a 2 - x 2max ))} l } 



F 2 i{gi,g 2 ,g 3 , [1 + exp(2a 2 (a 2 



Xlr, 



Here 5*1 is the S-matrix from the case IV ll and is given by the formula (|B10|) . The S-matrix related to the case IV 3 
(call it S 2 ) is related to Si through the condition (IB1[) . In the gluing formula (|B1|) S 2 stands for Ki and Si for Ki-i. 
Also 2 are the solutions of the Poeschl- Teller case 



Va(x) = e iux F 2 i [gi,g 2 ,g 3 , [l + exp(2a 2 (a; - x 2max ))] l ] , 



and 



9 a (x) =e- lux F 2 i [ 5l , 52 , 5 3,[l + exp(-2a 2 ( 



x - x 2r . 



0)] 



(B12) 
(B13) 



^(1-1)1.2 are the Morse potential solutions given as (|B2[) . 
We see that if there is a regime in which the modes do not 
feel the smaller peak, then they are given by the poles 
of Si and must be in the same time poles of S 2 . The 
formula for S 2 in this case (for frequencies that represent 
poles of Si) simplifies to 



S 2 = 



*(i-l)2*<l - (i-l)2 



i2*(i-l)2 



('-1)2 



(B14) 



I 

There is a possibility that for the case of the low modes 
in the two peak model one might be able to neglect the 
infinite depth of the valley between the infinite barrier 
and the Morse peak. In such case one can take a simpler 
model where the Morse potential directly follows after the 
infinite barrier. The Si term appearing in the formula 
(jBlip will be in such case given by a much more simple 
expression: 



Si (w,ai, V ) 



M{si,s 2 ,z(0)} 
MK, S2 ,z(0)}' 



(B15) 
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